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STATISTICAL ANALYSIS OF THE MIXED FRACTIONAL 
ORNSTEIN-UHLENBECK PROCESS 

P. CHIGANSKY AND M. KLEPTSYNA 


Abstract. Estimating the drift parameter of the Ornstein-Uhlenbeck process is a clas¬ 
sical problem in statistics. It can be approached in a number of ways, but only likelihood 
based estimators are known to attain the optimal accuracy in the large sample asymptotic 
regime. Construction and analysis of such estimators can be quite challenging, depend¬ 
ing on the properties of the driving noise. Essentially the only class for which a fairly 
complete picture is available consists of the fractional Brownian motions (fBm). In this 
paper we construct the Majcimum Likelihood Estimator (MLE), when the driving noise 
contains both standard and fractional components. Such mixtures arise in mathematical 
finance in the context of arbitrage modelling. We show how asymptotic analysis can be 
carried out in this case, using some recent results on the canonical representation and 
spectral structure of the mixed fBm. 


1. Introduction 


Estimating drift parameter 0 G M from a sample path of the Ornstein-Uhlenbeck type 
process = {Xt,t G [0,r]), generated by the equation 

Xt = Xo + 9 [ Xsds + Vt, t>0 (1.1) 

Jo 

is a prototypical problem in statistical inference of random processes. In its classical form, 
with the driving process V being the standard Brownian motion B = {Bt,t G [0, T]), 
this problem has been extensively studied since the late 60’s. In this case the probability 
measures /rj, 0 G M induced by are equivalent and the corresponding likelihood function 
is given by the Girsanov theorem: 


d^il 


(A^) = exp 


j-T 1 nl 

I / XtdXt - 

Jo ‘ ‘ 2 io 


X?dt 


( 1 . 2 ) 


Consequently, the MLE, being the unique maximizer of the likelihood, is given by the 
simple formula 


Ox = 


fo XtdXt 
/o" Xfdt ' 


T > 0. 


(1.3) 


It is known to attain the optimal asymptotic accuracy in the local minimax sense as T ^ oo 
and its limit behaviour is determined by the sign of the drift parameter 6. In the stable 
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case, corresponding to 0 < 0, the estimation error 6t — 0 is asymptotically normal at the 
usual parametric rate y/T : 

Vt(6t- e)N{0,26), V0<O (1.4) 


where the convergence is in distribution. Entirely different asymptotics emerges in the 
neutrally stable and unstable cases, 9 = 0 and 0 > 0. A comprehensive account of the 
early and more recent developments concerning the Brownian case can be found in [21], 
[19]. 

The drift estimation problem can be approached in a number of alternative ways, which 
typically produce reasonable estimators (see e.g. [14], [25]). However, without being based 
on the likelihood function, these estimators are asymptotically subefficient as T —>• oo, at 
best up to a finite gap with respect to the information bound. Construction and analysis 
of the likelihood based estimators, on the other hand, requires a convenient formula for 
the likelihood function, which can be hard to find for a given driving process V. 

The MLE for 9 can be constructed, taking advantage of the Girsanov theorem, if the 
process V is Gaussian and it admits a canonical innovation representation, [12]. More 
precisely, suppose that there exists a pair of deterministic kernels g{s, t) and g{s, t), so that 

Vt = f g{s,t)dMs and Mt = f g{s,t)dVs, t>0 (1-5) 

Jo Jo 

and M = {Mt,t > 0) is a continuous martingale with a strictly increasing quadratic 
variation (M) = t > 0). Here the stochastic integrals are defined in some reasonable 

sense, e.g., through approximation by simple functions. 

Then integrating the kernel g{s, t) with respect to A, we obtain a semimartingale: 

g{s, t)dXs = 9 [ QsiX)d{M)s + Mt 
Jo 

where 

The filtrations generated by X and Z coincide, = T^, t > 0, and by the Girsanov 
theorem the measures , 0 E M are equivalent with the likelihood function of the form 
(cf. (1.2)) 

||(A^) = exp (9 £Qt{X)dZt - £ Qt{Xfd{M)t^ . 



The MLE is therefore given by (cf. (1.3)): 

Jo 

The principle difficulty with this method is that the canonical representation (1.5) may 
not exist or can be hard to find in a suitable form for a given process V. Moreover, even 
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if such a representation is available, it may not be immediately clear how to approach the 
asymptotic analysis of the estimation error 


6t — d = 


Qt{X)dMt 


(1.7) 


Qt{Xyd{M)t 

Consequently likelihood based estimators have been studied only for a few processes beyond 
the standard Brownian motion. 

This program has been realized in [17] for the Ornstein-Uhlenbeck process (1.1), driven 
by the fractional Brownian motion ,t > 0), that is, a centered Gaussian process 

with covariance function 




1 


= Ut^^ + s^^-\t-sr), s,t>0, 


where H G (0,1] is the Hurst parameter. For H = ^ the fBm coincides with the standard 
Brownian motion, but otherwise has quite different properties, see e.g. [11], [22], [4]. 
In particular, for H > ^ it’s increments exhibit long range dependence, which makes 
fBm an important modelling tool (see, e.g., [3]). The MLE in [17] is shown to satisfy 
the asymptotics (1.4), using the canonical representation for the fBm, also known in the 
literature as Molchan-Golosov transformation (see also [23], [15]). 

A large class for which the canonical representation is available consists of mixtures 


Vt — Bt + Gt, t > 0 


( 1 . 8 ) 


where B is the standard Brownian motion and G is an independent Gaussian process. Such 
processes emerge in engineering applications, including signal detection problems (see, e.g., 
[16]). The kernels g{s,t) and g{s,t) in this case were constructed in [13], using solutions of 
certain integral equations, under the assumption that the covariance function of G satisfies 


EG.G. = 



K{u, v)dudv 


0 JO 


with a square integrable kernel 


fT f-T 

/ / K{u,v)‘^dudv < oo. 

Jo Jo 


(1.9) 


( 1 . 10 ) 


A natural question in this regard is whether the asymptotics (1.4) remains intact for such 
processes? 

In this paper we will answer this question affirmatively for the mixture with Gt '■= Bj^, 
H G (0,1), often referred in the literature as the mixed fBm. The interest in this process 
has been triggered by the paper [6], which revealed a number of its curious properties, 
useful in mathematical finance (see [7], [2]). Further related results can be found in [8], [1], 
[27], [5], [10]. 

We will show that the large sample behaviour of the MLE in (1.6) is governed in this 
case by a certain singularly perturbed integro-differnetial equation. While some elements 
of our approach remain applicable to more general mixtures of the form (1.8) (see Remark 
2.2 below), the results also indicate that the ultimate answer to the posed question requires 
quite delicate estimates on its solution, which might be hard to obtain in general. 
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2. The main result 

Consider the mixed fractional Ornstein-Uhlenbeck process (1.1) with V = B + , 

where B and B^, H E (0,1) \ {|} are independent standard and fractional Brownian 
motions. The conditions (1.9)-(1.10) are satisfied in this case only for > | and therefore 
the innovations construction from [13] cannot be used for the whole range of H. Instead 
we will use an alternative canonical representation suggested in [5] based on the martingale 

Mt:=E{Bt\3'Y). (2.1) 

To this end consider the integro-differential Wiener-Hopf type equation: 
d 

g{s^t) + — / g{r,t)H\s — sigii{s — r)dr = 1, 0<s^t<T. (2.2) 

ds Jq 

By Theorem 5.1 in [5] this equation admits a unique solution for any H E (0, ^]. It is 
continuous on [0,T], and the martingale, defined in (2.1), satisfies 

Mt= [ g{s,t)dVt and {M)t = f g{s,t)ds, fE[0,r], 

Jo Jo 

where the stochastic integral is defined for L^(0,T) deterministic integrands in the usual 
way. By Corollary 2.9 in [5], the representation (1.5) of V in terms of M holds with 

and consequently the MLE of 0 is given by (1.6). 


Theorem 2.1. The MLE for the drift parameter of the stable mixed fractional Ornstein- 
Uhlenbeck process satisfies (1.4). 


Remark 2.2. In the course of the proof we will show that the assertion of this theorem 
holds if the following two conditions are met: 



d , d 


2 

dt < oo. 


(2.4) 


and 


1 

— max 
t 


/ dt d{M)t \ t^oo 


AO. 


\d{M)fi dt J ' 

This part of the proof can be carried out in a greater generality, at least for sufficiently 
regular G with stationary increments. The equation (2.2) in this case takes the form 


9{s,t) + -^[ g{r,t)-^v{s-r)dr = l, 0<s^t<T, (2.6) 

as Jq dr 

where v{t) := EGfi Assuming that it has the unique solution, such that the aforementioned 
construction remains valid, the conditions (2.4) and (2.5) would still guarantee asymptotic 
normality (1.4). 

The second part of the proof amounts to verifying these conditions and this is where 
the particular structure of the fBm comes into play in an essential way. We show that the 
growth of the martingale bracket and its derivatives is governed by the solutions of the 
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singularly perturbed equation, obtained from (2.6) by an appropriate scaling. In the fBm 
case with v{t) = , asymptotic analysis of this equation is carried out using spectral 

approximations, recently obtained in [9]. It is not clear at this point, how similar results 
can be derived in the more general case. 


3. Proof of Theorem 2.1 


The proof is inspired by the approach in [18]. In view of (1.7), in order to prove (1.4) it 
is enough to show that (Theorem 1.19 in [19]) 

rT 


II 


We will check the stronger condition in terms of the Laplace transform 


£r(lt) := Eexp ( —fj, 


1 


Qtd{Mh 


T—>-oo 


exp 


m)' 


fJj G K+5 


(3.1) 


which in addition to the weak convergence (1.4), also implies the convergence of moments. 
The main difficulty in implementation of the method from [18], is that in the setup under 
consideration the kernels g{s,t) and g{s,t) do not admit explicit formulas. 


We will split the proof into several lemmas. Our first lemma shows that the process Q 
can be expressed as a stochastic integral with respect to Z, whose integrand involves the 
derivative d{M)t/dt. This derivative exists and is continuous by Theorem 2.4 in [5]. 


Lemma 3.1. 


where 



ii{s,t)dZs, 




1 / dt ds \ 

2 VWi ^ d{M)J ■ 


Proof. By Corollary 2.9 in [5], = Tf and 

Consequently, 



with g{s,t) given in (2.3). 


Qt — 


d 


d{M), 

d 

''d{M), 


[ 9{s,t)Xsds =[ g{s,t) [ g{r,s)dZrds 

Jo Jo 

!.([ s^ds^ dZ'p — ip {t ^ dZ^ 


with ^ 

where the equality f holds, since the integrand vanishes at r = t. Note that '0(s, t) does 
not depend on 6 and thus we can set 0 = 0 in all further calculations. For this choice 
Qt = Jq 'ip{r,t)dMr and hence 

KQtMs = f 'ip{r,t)d{M)r, s<t. 

Jo 










6 


P. CHIGANSKY AND M. KLEPTSYNA 


It follows that 


and 


s rt 


/O ^0 
r*s rt 


Define 


then 


si’-.t)Vrdr s(r,s)dVr, 

/ 9 {T,t)VrdT f g{r,s)dVr=Ef f g{T,t)g{r, s)^EVrVrdrdT = 

Jo Jo Jo Jo 

gir, s)g{T, t)^(TAr + ~ ~ '^\‘^^)dTdr = 

g{r, s)g{T, t) — H\r — rp'^“^sign(r — r)^ drdr = 

r s rt r s rt 

/ 9{r,s) / g{T,t)dTdr + H / g{r, s)r‘^^~^dr / g{T,t)dT- 
Jo Jr Jo Jo 

/ fi'(Ds) / ^(p, t)|?’— Tp^~^sign(r — r)(ir(ir. 

Jo Jo 

4>ir,t) ■= j 9 {T,t)dT - H J g{T,t)\r - sign{r - t) dr, 


10 JO 
r*s 


^ 4>{r,t) =-g{r,t) - H-^ [ g{T,t)\r - ^sign(r - r)(ir =-1 


dr 


dr 


and 


(f>{r,t) = (l){0,t) - r = f g{T,t)dT + H f g{T,t)T‘^^ ^dr - r =: {M)t + H{N)t - r. 
Jo Jo 


Gathering all parts together, we get 
, , d d 


d{M), d{M)t V 


(H{M)t{N)s+ / g{r,s)(^{M)t + H{N)t-ryr\ = 


jj d{N)s ^ d 


t + H{N)t 


d{M) 


g{r, s)dr — 


s Jo 


d dr 
d{M)s d{M)t Jo 


g{r, s)rdr = 


d{M), d{M)t 

„:M± 4. I 4. r, d{N), 

d{M). ^ d{M),‘ 

Integrating the equation (2.2) gives 

[ g{s,t)ds+ [ -^ [ g{r, t)H\s — r\‘^^~^ sigii{s — r)drds = t, 

Jo Jo ds Jo 

or, equivalently, {M)t + 2H{N)t = t, where we used the symmetry g{s,t) = g{t — s,t). 
Consequently, 

d{N)t _ dt 


1 + 2H 


d{M)t d{M)t' 


M. f) = + 1 + ff dWt ^ 1 / dt ds 

^ ^d{M), + ^d{M)t 2 [d{M)t ^ d{M), 


and 
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□ 

Lemma 3.2. Assume that {M)t is twice continuously differentiable. Then the convergence 
(3.1) holds under the conditions (2.4) and (2.5). 

Proof. By Lemma 3.1 

Qt= [ 'f{r,t)dZr = ]-ii{t,f)Zt + ]: [ 'f{r,r)dZr. 

Jo 2 I Jq 


Denote Y) = / 'f{r,r)dZr,ih.en 

Jo 

dZt = - \e\Qtd{M)t + dMt = -]^\e\i,{t,t)Ztd{M)t - ]^\e\Ytd{M)t + dMt = 




dWt, 


where W = iyVt,t G [0, T]) is a standard Brownian motion. Similarly, 


dYt = -\ 


1 , 


2, ^^(t, t)Ztd{M)t - \0\'if>{t, t)^Ytd{M)t + -fit, t)dMt = 

- ]^\6\'f{t,t)Ztdt - \9\]^Ytdt + ■s/Tp{t,t)dWt. 

The vector Q = {Zt,Yt)~^ solves the linear system of ltd stochastic differential equations 


dCt = -YA{t)Ctdt + b{t)dWt 


with 


1 


1 


A{t) = 

1 

The Laplace transform in (3.1) reads 

rT 


'if{t,t) I and b{t) = ( ^ /^{t,t) 

VWJ). 


Tr(l^) = Eexp Q^d{M)t'^ =Eexp J {'il;{t,t)Zt + YtYd{M)t 


t - 


Eexp 

where 


[Vi>{t:t)Zt + ^^ Yt'j dtj =Eexp^-^^iy C^R(t)(tdt^, 


m = 


''f{t,t) 
1 


V'(Li) 

By the Cameron-Martin type formula (see Section 4.1 of [18]), 

rT 


Lriti) = exp - 


4r 


f tr (r{s)R{s)^di 

Jo 


where r(t) is the solution of the Riccati equation 


|0| 


1^1 


r(t) = -JAA{t)T{t) - JAr{t)A{t )' + B{t) - ^T{t)R{s)T{t), t G [0,T] 


2T 


(3.2) 


(3.3) 
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with B{t) = b{t)b{t)~^ , subject to the initial condition r(0) = 0. 

The equation (3.3) has a unique positive definite solution for any E M and all T large 
enough, which can be expressed as the ratio T(t) = where and <h 2 (i) 

solve 

subject to <l>i(0) = I and <^ 2 ( 0 ) = 0. Since ti A(t) = 2 and 


(3.4) 


1^1 




it follows that 


$r^(t)$i(t) = Y^(t) + ^r(t)i?(f) 

^tr {r{t)R{t)) = tr (4>j"^(t)4>i(t)) - |6»|. 
By Liouville’s formula tr (<f>j"^(t)<hi(t)) = ^logdet (4>i(t)) and 


JL 

2T 


(3.5) 


tr {T{t)R{t))dt = logdet ($i(r)) - |0|T. 

Let us now calculate logdet (<hi(T)). To this end, define 

H ";) 

and note that R{t) = JA{t), B{t) = A{t)J and JA{t)J = A{t)~^ . If we set ^ 2 {t) ■= 
and multiply the second equation in (3.4) by J from the left, we obtain the system 


^,it)=^%,{t)A{t) + ^Ut)A{t) 
^2{t) = <^i{t)A{t) -^-^^2{t)A{t), 


(3.6) 


subject to ^i(O) = I and ^ 2 ( 0 ) = 0. The matrix 

2 , 

has two real eigenvalues ±Ar with At = (f)^ + ^ and the corresponding eigenvectors 


M JL 
2 27 

1 -ti 


V+ = I 


V_ = I 


where = -^ ± Diagonalizing (3.6) we get 

$i(t) = a;)iTi(t) + a;^T2(t), 


where 


ti(t)= XTTi{t)A{t) 

T2it) = -XTT2it)A{t) 


(3.7) 
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subject to Ti(0) = —T2(0) = \/\J (|)^ + Therefore 

logdet (<f>i(T)) = logdet (a;J;Ti(t) + a^T 2 (t)) = 

logdet + logdet + ^Tj'^(t)T 2 (t)^ . (3.8) 

Note that = |0| + 0{T~^) and detTi(O) = 9~‘^ + 0{T~^) and thus 

log det (ayTi(t)) — |0|r := log(o^)^ + log det Ti(0) + f At tr — |0|r = 

Jo 

2TXt - \e\T + 0(T-^) = T + ^ _ |0|^ + O(T-i) = ^ + 0{T-^). 

The desired limit (3.1) now holds by (3.2) and (3.5), if we check that the last term in (3.8) 
vanishes: 


|T^^(r)T2(r)|| 0, 


or, equivalently, 


T||T-i(T)T2(r)|| ^^0 


(3.9) 


since |a^/oJ| ~ T~^ as T —^ oo. Here and below we use the operator norm for matrices. 

We will need an estimate for the growth rate of the solution of the second equation from 
(3.7). To this end, denote gt := a-nd fix an arbitrary vector v E then 

vj ■= V^T2{t) 


9t 0 

0 l/gt 

satisfies the ODE vt = H(t)vt with the symmetric matrix 


H{t) = 


_ (-Xr + gt/gt 


—Xt 


—Xt 

-At - gt/gt) ’ 


The maximal eigenvalue of H{t) is 

-At + \/a| + {gt/gtf < ^ {gt/gtY 

and thus, under the assumption (2.4), 


\vt \\2 < Ibilb exp 


1 


(gt/gt) dt] < ||ui|| 2 C'i, r> 1 


.2At Ji 

with a constant Ci independent of T. Therefore 

||T2(r)|| < C2max(5rr, I/st), 

with a constant C 2 , which depends only on the initial condition T 2 ( 0 ). 
Further, note that 


dt 


Tr^(t) = -ATAi(t)Tr^(t) 
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which under transposition and mnltiplication by J from the right becomes 

^(Tr^(t)J) = -\TT^^{t)A^{t)J = -\TT^'^{t)JA{t)JJ = {t)J)A{t), 

i.e. J and T 2 (t) solve the same equation with possibly different initial conditions. 

Therefore 

||Tr'(T)T2(r)|| < ||Tr'(r)J||||T2(r)|| <C3max(5|.,l/g2), 
with a constant C^. The claim of the lemma now follows from (3.9). 

□ 

It is left to check conditions of Lemma 3.2, which we do separately for iL > ^ and 

H<1. 


Lemma 3.3. For H > 


— {M)t ~ T 


1 

2 

1-2H 


and 


d , d 

-r- log 

dT ^dT^ ' 


- 1—2 


as T —>■ oo 


and thus the conditions of Lemma 3.2 hold. 


Proof. For H > ^ the derivative and integration in (2.2) can be interchanged and it takes 
the form of integral equation 

g{s,t) + f g{r,t)cH\s — r\‘^^~‘^dr = I, 0<s<t<T, 

Jo 

where ch = H{2H — 1). By Theorem 2.4, [5] in this case 

{M)t= [ g‘^{t,t)dt, 

Jo 

with g(t, t) > 0 for all t > 0. 

Define small parameter e := A simple calculation shows that the function 

Us{x) := T‘^^~^g{xT,T) solves the integral equation 


/■ 


eue{x)+ / CH\y-x\^^ ^Ue{y)dy = l, a: e [0,1] 


and, moreover. 


d{M)q 

dT 


= g\T,T)=e^uUl)- 


(3.10) 


(3.11) 


For any e > 0 the equation of the second kind (3.10) has a unique solution, continuous 
on the closed interval [0,1] (see e.g. [26]). For e = 0 it degenerates to the equation of the 
first kind, whose unique solution is known in a closed form [20]: 


uo{x) = anx^ ^{l — x)^ 


-H 


X G (0,1) 


where an is an explicit constant. Note that uq explodes at the endpoints of the interval 
and therefore it is reasonable to snggest that Ueil) ^ oo as e ^ 0. 
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To estimate the growth of rie(l) we will use the following asymptotic approximations 
for the ordered sequence of eigenvalues and the scalar products with the corresponding 
eigenfunctions for the integral operator in (3.11) (see Theorem 2.3, [24]): for H > ^ 

\n ~ 

i^n) ^ n —oo (3.12) 

(1, ¥52n+l) ~ ^ 


where Cj’s are nonzero constants and {3 G (0,1). The eigenfunctions with even indices are 
antisymmetric around the midpoint of the interval and hence (1, </? 2 n) = 0. 

Taking scalar product of both sides of (3.10) and using these estimates we obtain: 


Us{l) = e ^ [ {uoix) - Ue{x))cHil - = 

Jo 

1 

/ J^nix)cH{l - xf^-‘^dx 
^n{e + K)Jo 

V 

Jo 


E 

n odd 




r-oo ^- 2 -H 


e + n 


h e + X 


dx ~ e 


oo y^-2 


+ 1 


dy, 


which yields the claimed asymptotics 


d{M), 

dT 


= e\2(l)~£ = r 


1-2H 


(3.13) 


The second condition is verified similarly: 

A. 

de 

-E 


^ /■i\ ^ i^nA) f \ \ 

''^£(1) = X ' N / Vn{x)CH\ 

n odd 


X - Ip'^ “^dx ~ 


-3;-H 


n 2 


n odd 

and hence 


(e + 


/ 


oo x-2-^ 


(e + 


dx ^ —s 2 


f 


,.3/1-2 


(y2H-l I'j 


2dy, 


^logT(M),)' = (^,„g,(T,T))E (^logT—„.(!))' < 


) =^ + 


rp2 


dT de 


dT 


1-2H A 
de 




dT Me(l) 


1 

2^2 ■ 


□ 


Lemma 3.4. For H < 


A 


{Mh 


const, and 


d , d 

^log^(A/)r 


T ^ as T ^ oo 


and thus the conditions of Lemma 3.2 hold. 
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Proof. For H < ^ the equation (2.2) takes the form (see Theorem 5.1 in 

A-2H 


where 


9 _ nzT 

CHmis) + = cnmis), s G (0,t], 


{'^g){s,t) = —2H— / g{r,t)r^~^ {r — s)^~2 dr 
J S 


Moreover, by Theorem 2.4 in [5], 


f{r)r 2 ^(s—r )2 ^ dr. 


i-H. 


ds 


{M)t = 


2-2H 






and it follows from (3.14) that 
p{t,t) = 


Ah ,H-i 
t ^CH 


2-2H 


CH 


]/ 2 i'^- 9 {pt))r 2 ^{t-r) 2 ^dr. 


(3.14) 


2-2H 

Let e := and define Ui;{u) := g{uT,T), u G [0,1], then 

f ^ 

p{T,T) = Ce ~2 / (1 — tte(x))u2“'^(l — x)“2“^(ix, 

Jo 

with a constant C > 0. The function Ue solves the equation 

eus + Kjj^Ue = (3.15) 

where Kh stands for the operator in (2.2). For H < ^ the inverse Kfj^ turns to be an 
integral operator with a certain weakly singular kernel (see (iv) of Theorem 5.1 in [5]). 
The limit equation is obviously solved by uq = 1 and hence 

/*1 

p{T,T) = Ce~2 / (^uo{x) — Us{x))x2~^(1 — x)~2~^dx. 

Jo 

Since uq, G L^(0, 1), 

Uo-U^= 2^(1, Pn)Pn - 2^ 


Pn 


' 'y Pn)Pn 2 ^ 


£ + An ^ 

An^(l,¥5n)_ _ e /T _ 

1 Pn — 2_^ i{^,Pn}Pn 

£ + An „ e + An 


Define h{u) := Cu 2 ^(1 — u) 2 ^^ then 

p{T,T) = e^'^ - ^{^,Pn){h,Pn)- (3.16) 

„ e + An 

By Theorem 2.3, [24], for iL < ^ the eigenvalues satisfy the same asymptotics as in (3.12) 
and therefore form an increasing sequence, in agreement with the fact that in this case 
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the operator Kh is not compact. Also we have ((•) ^h,ipn) rsj Chn^ ^ for any j3 G (0,1) 
and any continues function with /i(0) 7 ^ 0. However the asymptotics of the averages for 
symmetric eigenfunctions is completely different in this case: {l,ipn) ~ C' 4 n“^. 

Now we can estimate the growth rate of the series from (3.16): 


1 


r{e)= ^ ~ 

£: “T A' 


-^+H 


n 2 


e + n- 


2H-1 


~ £ 2 


00 y ^ H 


n odd n odd 

and hence 

p{T,T) ~ const., T —)• 00 , 

which verifies the condition (2.5). 

Further, differentiating the series in (3.16), we get 


Jo 


- 2/1 + 1 


dy 




ds ^ £ “ 1 “ A' 

n 

E 


n 2 +^ 


(e + n2//-i)2 


~ —E 2 


yk-'iii 


r y^- 

Jo 


+ iy 


rdy, 


and 


— log—(M)^) = 
dT ^dT^ ' ' 


A 


logp(r, T) 


fded 1 


T 


-2 


which verifies (2.4). 

Remark 3.5. In the simpler regression problem 

Xt = et + Vt, te[o,T] 


□ 


the MLE of 0 G M is given by 

Consequently the estimation error is normal with zero mean and its variance is controlled 
by the growth rate of {M)t, rather than the derivative d{M)T/dT as in the Ornstein- 
Uhlenbeck problem. Finding asymptotics of the bracket {M)^ amounts to singular per¬ 
turbation analysis of the equations (3.10) and (3.15) with respect to weak convergence (cf. 
(3.13)), which can be carried out either directly (see the discussion concluding Section 7.1 
in [24]) or using the spectral asymptotics as above. 

The corresponding limit variance is 


E(0 t - of ~ 


vhT^^-^ 


H> ^ 

H<1 


with vh 


2HT{H + - 2H) 

r(i-^) 


and it follows that the asymptotic is dominated by the fractional component for 77 > ^ 
and by the standard Brownian component for 77 < + 
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